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ABSTRACT 

We present a simple semi-analytical model of nonlinear, mean-field galactic dynamos and use 
it to study the effects of various magnetic helicity fluxes. The dynamo equations are reduced 
using the 'no-z' approximation to a nonlinear system of ordinary differential equations in 
time; we demonstrate that the model reproduces accurately earlier results, including those 
where nonlinear behaviour is driven by a magnetic helicity flux. We discuss the implications 
and interplay of two types of magnetic helicity flux, one produced by advection (e.g., due to 
the galactic fountain or win d) and the other, arising from anisotropy of turbulence as suggested 
by IVishniac & Choi dMTh . We argue that the latter is significant if the galactic differential 
rotation is strong enough: in our model, for R^j < — 10 in terms of the corresponding turbulent 
magnetic Reynolds number. We confirm that the intensity of gas outflow from the galactic 
disc optimal for the dynamo action is close to that expected for normal spiral galaxies. The 
steady-state strength of the large-scale magnetic field supported by the helicity advection is 
still weaker than that corresponding to equipartition with the turbulent energy. However, the 
Vishniac-Cho helicity flux can boost magnetic field further to achieve energy equipartition 
with turbulence. For stronger outflows that may occur in starburst galaxies, the Vishniac-Cho 
flux can be essential for the dynamo action. However, this mechanism requires a large-scale 
magnetic field of at least ~ 1 to be launched, so that it has to be preceded by a conventional 
dynamo assisted by the advection of magnetic helicity by the fountain or wind. 
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1 INTRODUCTION 

' Conservation of magnetic helicity has recently been recognized as 
\ a key constraint on the evolution of cosmic magnetic fields, es- 
. pecially those produced by the large-scale dynamo action. Any 
' large-scale magnetic field (if not supported by external currents) 
must have both poloidal and toroidal parts in order to mitigate 
diffusive losses; therefore, any self-sustained large-scale magnetic 
field must have non-zero magnetic helicity. One of consequences 
of the helicity conservation is the suppression of the a-effect in 
a medium of high electrical conductivity. Since the large-scale 
magnetic field has a non-zero helicity in each hemisphere through 
the mutual linkage of poloidal and toroidal fields, the dynamo 
must be producing small-scale helical magnetic fields with the 
opp osite sign of the magnetic helicity. In several closure mod- 
els ( Pouquet, Frisch & Leorat 1976; Kleeorin &Ru zmaiki3l 19821 : 
iGruzinov & Diamond 1994 ; Blackman & Field 200i ). this leads to 
a suppression of the dynamo action due to the current helicity of the 
growing small-scale magnetic field. The suppression can be catas- 
trophic in the sense that the steady-state energy density of the mean 
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magnetic field is as small as 7?~^ times the kinetic energy density, 
where Rm S> 1 is the magnetic Reynolds number. Thus, efficient 
action of the mean-field dynamo requires that the magnetic helicity 
due to the small-scale magnetic field is removed from the dynamo 
active region (Blackman & Field.2000. . .2001. : . Kleeorin et alj200oh . 



Several mechanisms have been suggested to produce 
the required helicity flux, including the ani sotropy of the 
turbu l ence produced by differential ro tation jvishniac & Choi 
l200ll : ISubramaman & Brandenburg|l2004l) , and t he non-uniformity 
of the Q-effec t ('Kleeori n et alj |2000|. |2002|) . More recenfly, 
IShukurov et al. suggested a simpler mechanism involving 

advection of small-scale magnetic fields (together with their helic- 
ity) by an outflow from the dynamo-active region, e.g., the galactic 
fountain or wind in the case of spiral galaxies. We explore these 
effects using a simple model where dynamo equations are reduced 
to a system of ordinary differential equations, which captures the 
salient features of the nonlinear dynamo action and helicity evolu- 
tion. The simplicity of the model allows us to gain deeper insight 
into the role and interaction of various mechanisms, and to more 
extensively explore the parameter space. We consider both an ad- 
vective helicity flux and that arising from the combined action of 
velocity shear due to differential rotation and anisotropy of the tur- 
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bulence produced by it. We demonstrate the efficiency of the galac- 
tic dynamo action provided it is facilitated by the helicity fluxes. 



2 HELICITY BALANCE AND MEAN-FIELD DYNAMOS 

The generation of a large-scale magnetic field by small-scale 
random motions is desc ribed by the mean-field electrodynamics 
( lKrause&Radleilll98d) . where relevant physical quantities are 
split into mean and fluctuating parts, e.g., B = B + b for mag- 
netic field and A = A ~\- a for the vector potential, where overbar 
denotes relevant averaging and b = 0, a = 0. This results in the 
mean-field dynamo equation 



dB 
~dt 



= \7 X {U X B + £ -rjj) , 



(1) 



where U is the mean velocity field, 77 is the Ohmic magnetic diffu- 
sivity, and J = \/ x B/fio (for the vacuum magnetic permeability 
we assume /io = 1 hereafter). Furthermore, 



£ = u X b — aB — TjtJ 

is the mean turbulent electromotive force (emf) (assuming 
isotropic turbulence), with a and rjt the turbulent transport 
coefficients responsible for the a-effect and turbulent magnetic 
diffusion, respectively, and u is the small-scale velocity field. The 
dependence of a on B and b is a subject of ongoing research, 
and our knowledge of this is based on limited closure models and 
numerical simulations. In s ome cl osure models, such as EDQNM 
( Pouguet. Frisch & Leorall Il976| ) and the r- approximation 
tBlackman & Field I2OO2I: iRMler. Kleeorin & Rogachevskii 
I2OO3I : [Brandenburg & SubramaniaiJ l2005bh . the effect of the 
small-scale magnetic field on the a-effect is described by 
terandenburg & Subramaniaij|2005ah 

Q = qk + am , (2) 

where ok represents the kinetic Q-effect, ok ~ ■ V x u, 

and Om = ^P~^Tj ■ b is the magnetic contribution to the a-effect, 
with p the fluid density and r the correlation time of the turbu- 
lent velocity field u (assumed to be short). Further, rjt — l^u^. 
It is not clear how general is the form Using a low Reynolds 
number approximation, where u nambiguous analysis is feasible, 

Subramanian & Brandenburg! ( |2007[) argue that Eq. lO is ac- 
ceptable, with aK modestly affected by magnetic field. 

To constrain a^i, we use the helicity conservation equation 
written in terms of the helicity density x of the small-scale mag- 
netic field, defined in a gauge-invariant form in terms of the num- 
ber density of the links of b. Since the small-scale magnetic field 
has finite correlation length, one can meaningfully introduce its 
helicity density and then derive its transport equation in the form 
jSubramanian & Brandenburg||2006l) 



1 + '- 



2£ B 2r]j ■ b, 



(3) 



where -F is the helicity flux density (specified below) and j = V x 
b. For practical purposes, x is approximately equal to a • 6 given 
that a is defined in the Coulomb gauge. 

In the steady state, Eq. {3} yields £ ■ B — — ■ F — r]j ■ b. 
For = 0, it follows that f-B^Oasr;— >0 under reasonable 
assumptions about the spectrum of current helicity j ■ b (specifi- 
cally that it peaks at a scale independent of 77 or at least at such a 
scale that r;j • 6 — > as 77 ^ 0). Hence the component of the emf 
parallel to the mean magnetic field vanishes, and so the dynamo 



becomes inefficient. This catastrophic quenching of the a-effect in 
highly conducting medium is, however, avoided if f 7^ 0. The 
simplest contribution to the flux. 



(4) 



arises from the net effect of advection of magnetic fields by a ve- 
locity fiel d U directed away from the dynamo active region, as sug- 
gested bv lShukurov et al. 1 2006h. We als o include another helicity 
flux proposed bv lVishniac & CholfeOOlh . 

In order to relate to x, we argue that the main contribution 
to Om comes from the integral scale of the turbulence lo = 2n/ ko 
dShukurov e7Zll2006l) , so that j ■ b ~ Ig^a ■ b and then 



3 ;2 
'0 



(5) 



To justify this relation, we note that the spectral density of am oc 
rj ■ b differs from that of x — a • 6 by a factor r^k'^. When 
the spectral density of x decreases with k fast enough, the spec- 
trum of am decreases with k as well, and then both x ™d 
Om are dominated by the integral scale, 2n/ko, which is the 
meaning of Eq. ([5j. For example, for the Kolmogorov spectrum 
Mk = k~^bl oc k~^^^, we have bk oc fc"^/^ and Tk oc 
k~^'''^, so that the spectrum of x is dfebfe oc k~^''^, that of 
the current helicity is jkbk oc k^^'^, and Tkjkbk oc k^^''^ for 
am. Then both am and x ^6 dominated by the integral scale. 
These arguments are similar to th o se pre sented in Sec t. lO.IV of 
'Zeld ovich, Ruzmaikin & So kolof^ fm^. Analytical jBlackmanl 
2003) and numerical results (Brandenburg & SubramanianI 2005bh 
suggest that the spectrum of the current helicity in fact decreases as 
k^^^'^; this lends further support to Eq. 

Introducing a reference magnetic field strength Bfq = pu'^ 
and defining the magnetic Reynolds number as Rm ~ rjt/V' we 
rewrite Eq. ((S}, using Eqs (HJl and ([S}, as 



dt 



27?, 



£ B Oir 

;2 I d2 ' n 

\ ^cq Jl-rr 



V- am?7 



(6) 



This equation should supplement the mean-field dynamo equa- 
tions to provide our model of nonlinear dynamo action. 
Since galactic discs are thin, it suffices to consider a one- 
dimensional model, retaining only the z-derivatives of the vari- 
ables (Ruzmaikin, Shukurov & Sokoloff 1988). In terms of cylin- 
drical coordinates [r, tp, z), we consider a mean flow consisting of 
differential rotation and vertical advection, with U — {0,U z)- 
Altogether, the system of nonlinear mean-field dynamo equa- 
tions has the dimensionless form 



dt dz 



dz2 



^ =R^Br-Ru^iU.BA + %^ 
at oz ^ ' oz-^ 



(7) 
(8) 



dc 



= -C (aB^ - R-^J ■ B + -Rv^ (amf/.) , (9) 



dt 



where, retaining only the ^-derivatives. 



dBr 
dz 



B, 



dB^ 
dz 



(10) 



and we have neglected the term proportional to Ra in Eq. ((8) thus 
adopting the acij-dynamo approximation for convenience; this term 
can easily be restored. Here the time and length units are /rjt and 
h, respectively, with h the semi-thickness of the disc; the unit of a 
is denoted ao, and that of U z is Uo; magnetic field is measured 
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Ru = 



^, R.^^, R.^^, C^2(^) 



(11) 



in the units of Bcq- We have introduced dimensionless parameters 
defined as 

Uoh 
Vt 

where G = rdD./dr and we have neglected Bz in comparison 
with the other two compone nts of the mean magne tic field, which 
is appropriate in a thin disc JRuzmaikin et alJI 19880 . We shall also 
use the dynamo number defined as D = RaRui- 



3 REDUCED DYNAMO EQUATIONS 

We reduce the dynamo equations to a dynamical system, where all 
the variables are functions of time alone, by using the 'no-z' ap- 
proximation ('Subr amanian and Mestellll993l . IMoss|[T995| . IPhillipsI 
[200 1). Specifically, the z-derivatives of all the variables are re- 
placed by appropriate division by the disc semi-thickness, d/dz 
1/h and /dz^ — » The field components Br, B^ appear- 

ing in the resulting ordinary differential equations can be thought 
of as representing either the corresponding mid-plane values or ver- 
tical averages. By its nature, the no-2; approximation captures the 
quadrupole modes of the dynamo which are dominant in a thin disc. 

The current helicity density of the large-scale field J ■ B van- 
ishes in the no-z approximation - see Eq. l llOt . Therefore, we calcu- 
lated it using a perturbation solution of one-dimensional, kinematic 
mean-field dynamo equations as described in Appendix IaI 



J B 



\(~-KDf'''BrB^ 



(12) 



When applying this estimate to study nonlinear sol utions, we recal l 
that D IX a = q:k + Qm. Furthermore, we follow Philligsl l l200lh 
who suggested that the accuracy of the no-z approximation can be 
improved by using 



dz2 



TT 



1 

ft2 



(13) 



These rescalings are consistent with the z-dependence of B ob- 
tained in Appendix |A] Derivatives similar to d{aB^)/dz have to 
be treated cautiously in the no-z approximation since this term is 
positive near z — and negative near z — 1 (note that a = at 
z = and B^ — at z — 1). The main contribution to the dynamo 
action comes from the vicinity of z = where jS^I is maximum, 
so that the par t of the above t erm important for the dynamo action 
is B^da/dz ( iSokolofllll996l) . where da/dz > near 2; = in 
galactic discs. This suggests that the no-2 approximation of this 
term is aB^/h rather than —aB^/h. 

The resulting reduced Eqs l|7](-([9]( are given by 



dBr 
dt 

dB^ 
dt 

dajn 
dt 



Ra (1 + Ctn 

TV 



\B4,-[Ru 



TT 

T 



Br 



— RuOtm — C 



Ru + —]B^, 



{1 + am){Bl + bI) 

Nl/2 



(i+am) -5 BrB^ + — 

O-ttct rin 



(14) 



(15) 



(16) 



where aK = ao is absorbed into Ra, so that a = 1 + am in terms 
of dimensionless variables, and Uz can now be put equal to unity. 
Note that D must be replaced by Z)(l + Qm) when Eq. l ll2b is used 
in the helicity evolution equation because D in that equation has 
the meaning of the instantaneous dynamo number. 




Figure 1. The phase portrait of the dynamical system jl4t -ll7t in the 
— 2 

{B , — am)-plane with and without vertical advection, ^ Rjj ^ 0.5. 
Parameter values used here ai'e Ra = 1, Ru: = —15, C = 50, 
Rui = 10^, the values of Ru are shown in the legend, and the time span is 

t ^ loo/iV^jt. 



Parameter values typical of the Solar neighbourhood of the 
Milky Way are as follows. The scale height of the ionized gas 
layer is ft ~ 500 pc, the integral scale of interstellar turbulence 
is ^0 =i 100 pc, which yields C = 2{h/lof ~ 50. With the turbu- 
lent magnetic diffusiv ity r;t — lO'^^cm'^s"^ and Uo ~ lkms~^ 
jShukurov et al.l l2006), we obtain Ru ~ 1.5; we consider a range 
Ru = 0-2. Furthermore, we use standard values of the dynamo 
parameters, Ra — 1 and R^: — —15. These estimates are often 
adopted as typical of spiral galaxies in general. However, many 
galaxies rotate faster than the Milky Way, and even within the 
Milky Way the angular velocity of rotation rapidly increases to- 
wards the centre. Likewise, the other parameters vary broadly be- 
tween galaxies and across a given galaxy. Hence, any identifica- 
tion of the above parameter ranges as generally applicable to spiral 
galaxies should be treated with caution. 

To illustrate the accuracy of the no-z approximation, we cal- 
culate the critical dynamo number (corresponding to dB /dt — 0) 
from Eqs J14t and il5l with a„^ = and Ru = as Dc ~ 
— (7r/2)^ ~ —9.6, as compared to the more accurate values (ob- 
tained from a numerical solution of the boundary value prob- 
lem) Dc ~ —8 and — 11 for qk = sinTrz and z, respectively 
( iRuzmaikin et al.|[l988[) . In the next section we demonstrate that 
relevant features of the n umerical solutions of t he nonlinear equa- 
tions 0-119), obtained bv lShukurov et al.l 1 I2OO6I) , are accurately re- 
produced as well. 



4 THE EFFECTS OF THE ADVECTIVE FLUX 

We solved numerically Eqs Ill4t-il6t firstly in order to verify the 
accuracy of the no-z approximation, and then to study the nonlinear 
evolution of the dynamo. The initial conditions used are 



Br=0, 



Bi 



10" 



-10" 



at f = 0. 



(17) 



The choice am|i=o = does not affect the results. 

In Fig. [T] we show the phase portrait of the system in the 
[B , —am) -plane, obtained for various values of Ru- For weak 
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6C 





I steady state 



Figure 2. The evolution of the magnetic field strength for various mag- 
nitudes of the vertical advective flux shows that magnetic field eventually 
decays for both small and large values of Rjj. Parameter values are the 
same as in Fig.[T] Time is given in the units of the magnetic diffusion time 
across the gas layer, about 7.5 X 10* yr. 



advection (Ru < 0.01), — Ofm monotonically increases with time 
to reduce the net Q-effect, so that the dynamo eventually becomes 
subcritical and magnetic energy reduces to negligible values. The 
dynamo action can resume again as soon as will have decayed 
together with the large-scale magnetic field and a = qk + Qm 
becomes supercritical again. Without the advective flux of helicity, 
the time scale for such a recovery of the dynamo would be of the or- 
der of the Ohmic diffusion time, as controlled by the term Om/Rn-i 
on the right-hand side of Eq. l ll6t . However, if Ru is of order unity, 
the relaxation time of the magnetic helicity becomes significantly 
shorter because of the advection represented by the first term on 
the right-hand side of Eq. J16b . and magnetic field grows again to 
exhibit nonlinear oscillations before it approaches a steady state. 
The time scale of the oscillations is long at about 2 x lO^'' yr for 
Ru = 0.05. The magnitude of B in the steady state first grows 
with 7?;;/, but then reaches a maximum at Ru ~ 0.3. Stronger ad- 
vection affects the dynamo adversely by removing B too fast. 

The variation of magnetic field strength with time is also 
shown in Fig.|2] The initial exponential growth of the magnetic field 
is catastrophically quenched in the absence of the advective flux, 
and the field decays at about the same rate as it grew (long-dashed). 
However, even a moderate advective flux (Ru ~ 0.3, solid) com- 
pensates the catastrophic quenching allowing the magnetic field to 
reach a steady-state value of about 0.1-Bcq. For stronger advection 
(Ru ^ 0.8, see the dash-doted curve), the dynamo action is sup- 
pressed again since the mean-field is removed too rapidly from the 
dynamo-active region. 

The above results ag ree quite accurately with those obtained 
by 'Shuk urov et alj j2006b who solved the differential equations 
([7]l-([9]l. This confirms the applicability of the no-z approximation 
to nonlinear dynamo equations. In addition, this indicates that our 
estimate of J • B in Eq. l ll2t is appropriate. 

To clarify further the effect of the advection on the dynamo, 
we consider the steady-state solution, d/dt — in Eqs l ll4l (-lll6b. 
where we introduce the critical dynamo number Dc such that, in 
the steady state, 



Then Eqs O and l[T5) yield 



D. = --(Ru + - 



Br 



\ n ) R^ 



(18) 



(19) 



the first of which shows that the critical dynamo number is affected 
by the advection; in this sense, the outflow hinders the dynamo ac- 
tion, as might be expected.The second of this relations shows that 
\Br\ <C 1-6,^1 if \Ruj\ ^ 1. ThenEq. il6i yields, neglecting B^. in 
comparison with B^, 



B 



i 

c 



where 



IV2 



(20) 



(21) 



— 1/2 

It is clear that B oc Rm for Ru ~ 0, which is the case of 
catastrophic a-quenching ( with B p erhaps exhibiting long-term 
variations around this level - iBranden burg & Subramanian 2005d). 
However, the magnetic field strength is approximately proportional 
to RU^ for Rm 2> 1 and small Ru- On the other hand, \Dc\ in- 
creases with Ru, so that there is an optimal value of Ru provid- 
ing maximum field strength. For D = —15, the steady state is 
nontrivial (B > 0) for Ru < 0.6, with B being maximum for 
Ru ~ 0.3. In the optimal case, Ru = 0.3, we obtain B ~ 0.1, 
where we recall that the unit magnetic field corresponds to equipar- 
tition between magnetic and turbulen t energy densities. We note 
that Eq. (HoJ is similar to Eq. (10) of IShukurov et all (l2006h . but 
here we use an arguably better estimate of J • B. 

The steady-state magnetic field remains of order C^^^^Bcq — 
0.1-Bcq, with C = 2{h/lo)'^ ~ 50, even in the presence of the 
advective flux. Therefore, we consider additional helicity fluxes to 
examine if they can lead to stronger magnetic fields such that B ~ 

Bca- 



5 THE VISHNIAC-CHO FLUX 



A flu x of magnetic helicity discovered by Ivishniac & Choi 
( I2OOII) relies on the anisotropy of turbulence which is nat- 
urally produced by velocity shear due to differential rotation 
('Sub ramanian & Brandenb urg|l2004l : IBrandenburg & SubramanianI 
20053). Indeed, the regular velocity shear Uy{x) produces an ad- 
ditional {/-component of the turbulent velocity via duy/dt ~ [u ■ 
\7)Uy, SO that the resulting anisotropy, produced during one corre- 
lation time T,is Uy/ux — 1 + rdUy/dx. The x and y directions 
can be identified with the radial and azimuthal ones in the galactic 
disc, so that dUy/dx is replaced by G = r dSl / dr. The anisotropic 
part of the turbulent velocity correlation tensor can be estimated 
as UxUy ~ tGu^, where u is the background isotropic turbu- 
lent velocity field. A convenient ex pression for the flux of Vish- 
niac & Cho has been obtained by Brandenbu rg & SubramanianI 
tOOSd) andlSubramanian & Brandenburg t200d) . Using Eq. (12) 
of IBrandenburg & SubrMnanianI i2005cl). with their Cyc c alcu- 
lated using the results of I subramanian & Brandenbur j ^20061) . we 
can represent the vertical flux of the magnetic helicity of the small- 
scale magnetic field in the following dimensional form: 



^iurfG{B 



^2 -=;2 



B, 



(22) 
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This flux win add tlie following term to the right-hand side of 
Eq. J16I ) written in the dimensionless form (and using the no-z ap- 
proximation): 



dam 
dt 



OF, 

dz 



dz 



(23) 



where dots denote the terms already included in Eq. ( I16t . Concern- 
ing the steady state, Eqs. dlSt and l |19t still apply and the steady- 
state strength of the mean field is given by Eq. ( I20t . but now with 



1 



1- |V2 + 7?£/(CDc) 



(24) 



instead of l l21t . Since < 0, the additional flux results in a 
stronger steady-state magnetic field than that with advective flux 
alone - cf. Eq. ( 121b . 

If \Rui \ is large enough, the additional term can formally lead 
to a singularity in B where the denominator in Eq. J24b vanishes 

2 

(and B < for larger values of \Ruj\). The reason is that the 
Vishniac-Cho flux produces > enhancing the hydrodynamic 
Q-effect. If i?J is large enough in comparison with C|-Dc|, the flux 
of Vishniac-Cho can lead to a dynamo action of its own, indepen- 
dently of the kinetic a-effect (Vishniac & Cho 2001). As shown in 
Sect. 15. ll equations governing this regime are in fact linear in B, 
hence B grows exponentially as t ^ oo, which corresponds to the 
singularity in Eqs l |20t and l l24t . Despite the formally linear gov- 
erning equations, this type of dynamo is essentially nonlinear and 
relies on either a strong seed magnetic field or an earlier conven- 
tional dynamo action to produce a strong magnetic field required to 
build up Qm at a sufficiently short time scale (Sect. lSTt . 

For the parameter values typical of the Solar neighbourhood, 
C = 50 and Dc = —10, the effect of the additional helic- 
ity flux becomes significant as R^^ approaches about —15. For 
\Rui\ 15, this effect does not lead to independent dynamo ac- 
tion (see Sect. l5.lt but rather increases the steady-state strength of 
the large-scale magnetic field to values close to Bcq. 

With the addition of the Vishniac-Cho flux, we plot in Fig. [3] 
the variation of the magnetic field strength driven by both helicity 
fluxes for a particular value of Ra and varying R^,, as obtained 
from a numerical solution of Eqs (I14b-J16b with the modification 
i23\ . The useful contribution of the Vishniac-Cho flux at larger val- 
ues of Ru is now evident. With the advective helicity flux alone, 
magnetic field decays if Ru is large enough; this is true with the 
Vishniac-Cho flux added, but now for larger values of Ru- For ex- 
ample, B decays for Ru > 0.8 if i?^; = —15 (either with the flux 
of Vishniac & Cho - Fig. [3^, or without it - Fig. O, but only for 
Ru > 1 if Rui — —17 (Fig. [3};). This is a result of the joint ac- 
tion of the two mechanisms, the QCLi-dynamo and that driven by the 
Vishniac-Cho flux: both become stronger as \Rui\ increases. For 
Rui ^ —15.5, the solutions for Ru = do not decay because 
the Vishniac-Cho flux drives a dynamo of its own (the dash-dotted 
curves in Fig.[3}5,c). 

Thus, the advective helicity flux facilitates the dynamo action 
by alleviating the magnetic helicity conservation constraint. On the 
other hand, it ensures that the dynamo action is eventually satu- 
rated. On the contrary, the helicity flux of Vishniac & Cho, also 
mitigating the helicity constraint, can lead to an unbounded growth 
of magnetic field and must be quenched by further physical pro- 
cesses. The advective flux can ensure such a quenching for moder- 
ate values of \Rui \ (e.g., for Ru = 0.3-0.8 in Fig.[3ll- 




Figure 3. Evolution of the magnetic tield strength with both the advective 
helicity flux and that of Vishniac & Cho, flux for various magnitudes of the 
velocity shear: (a)H„ = -15, (b) _R^^ = -16 and (c) i?„ = -17. The 
values of Ru presented are: (dash-dotted), 0.3 (solid), 0.5 (dashed), 0.8 
(long-dashed), 1.0 (dotted) and 3.0 (short-long-dashed). The Vishniac-Cho 
helicity flux supports the dynamo action for larger values of Ru '■ note decay 
at Ru = 0.8 in (a) but not in (c). The dynamo action is similarly facilitated 
for small Ru , where solution with Ru = decays in (a) but not in (b) and 
(c). The parameters used here are C = 50, Rm = 10^, B</)|t=o = 10~^ 
and Br\t—o = 0, with the values oi Ra and R^ indicated in the legends. 



5.1 Dynamo action due to tlie Vishniac-Cho flux 

In this section we present a n illustrative model of dynamo action 
driven by the helicity flux o f lVishniac&Ch ^faOOl'i alone. For this 
purpose, we use Eq. l|9) with additional term l l23t but now we do not 
employ the no-z approximation, take Ru ~ and neglect any hy- 
drodynamic a-effect, qk = 0. To make the system easily tractable, 
we make two key assumptions which can be justified a posteriori 

2 2 

using the solution obtained: we neglect B,. in comparison with B^ 
in Eq. Q and assume that J • i? is negligible. Neglecting also 
Qm/i?m, this yields 



dt 



iCBj. 



dz 



where the dimensional form of Fz is given in Eq. l l22t . and we 
use dimensionless variables as defined in Sect.|2] Then Eqs l[7]l-([9]l 
reduce to 



dBr_ 
~dt 

dB^ 
dt 

dOm 
dt 



d — 



+ ■ 



d^Br 

dz'^ 



RuiBr H 77-^ , 

az^ 

2 dB^ 

-amCBs + R 

OZ 



(25) 
(26) 
(27) 



where we have put Ra — 1, which is equivalent to choosing the 
unit of to be r]t/h. The characteristic time of the variation 
of Qfm, of the order of the turbulent time scale, is much shorter 
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than that of the magnetic field (see Sect.[6ll. Therefore, we can take 
dctn^/dt — 0, which yields 



.Ru. dlnB 



" c dz ■ 

Perhaps unexpectedly, Eqs. ( 125 1 and i26i now become linear. 



dt 
dt 



= -2- 



C dz^ 

RujBr 



dz^ 



dz^ 



(28) 



(29) 



and have the following exact quadrupolar solution satisfying the 
vacuum boundary conditions l lA3t 



04, 



1 



cos{ttkz) exp {"/t) , 



where growing magnetic fields have 



i + n , n = 0, 1, 



and K is an arbitrary constant. Given that Ru, < in galactic discs, 
the dynamo produces a growing magnetic field if R^^ < R^^ with 



(30) 



As we argue below, this is a necessary but not sufficient condition 
for the dynamo action of this type. 

The assumptions made to derive this solution can now be veri- 
fied: \Br/B4,\ = n/V2C < Iforthelowestmode, and J-B = 
since Br/B^ is independent of z in this approximation. 

An essential feature of this mechanism is that it cannot be 
launched unless the large-scale magnetic field is strong enough: 



for B ^ 1, we have 



<^ 1 and the field must decay. For 



am to grow, its variation rate CB should be larger than the mag- 
netic diffusion rate across the disc, 7r^/4 of Eq. il3l . This yields 
B > Tv/{2\/C) ~ 0.2, where we recall that magnetic field is mea- 
sured in the units of Bcq. For galactic parameters, the minimum 
magnetic field is of order 1 /iG. The initial magnetic field needs to 
be even stronger if we take into account that it first decays until am 
has grown enoug h. These arguments equall y apply to the helicity 
flux discussed bv lKleeorin etal] ( l200(ll2002h . 

Therefore, conditions for the dynamo action driven by the 
Vishniac-Cho alone are, firstly, the inequality l |30| l and, secondly, 
the initial large-scale magnetic field must be strong enough, say 
B |t=o > Bq, where Bo strongly depends on R^j. As a result, the 
dynamo threshold Ri^^ is a function of the initial magnetic field; in 
this sense, Eq. J30t is just a necessary condition. Numerical so- 
lution of the dynamo equations in the no-2: approximation with 
qk = shows that magnetic field decays for Ruj > —22 for any 
initial magnetic field and grows for R^j < —23 if Bo > 0.4 and 
for Ruj < —40 if Bo > 0.1. This dynamo mechanism seems to be 
suitable for an additional amplification of the large-scale magnetic 
field produced by the conventional mean-field dynamo (assisted by 
the advective helicity flux as described in Sect.|4l( closer to equipar- 
tition with turbulent energy. 

The Vishniac-Cho dynamo can be saturated by quenching the 

corresponding helicity flux, with dF^/dz in Eq. ( 1231 ) multiplied by 

1/(1 + B ). We have confirmed, using the no-z approximation, 

that this indeed produces a steady state with B = Oil). Curiously, 

2 

the 'standard' Q-quenching, with a^i am/(l + B ) in both 
Eqs ( |25t and ( |27| l, cannot lead to a saturated state because of the 
trivial cancellation of the quenching factor. 



6 CONCLUSIONS AND DISCUSSION 

The simple model suggested above reproduces a feature of the 
mean-field dynamo which has become well known: if the helicity of 
the small-scale magnetic field cannot escape from the dynamo ac- 
tive region, the mean magnetic field eventually decays to negligible 
levels being constrained by the conservation of magnetic helicity. 
Correspondingly, Eq. ^ yields B ~ R^^''^ < 1 for Ru = 0. 
An outflow from the dynamo region can prevent this catastrophic 
suppression of the dynamo as it carries away small-scale mag- 
netic fields together with their contribution to the total magnetic 
helicity, giving breathing space to the large-scale magnetic field 
dShukurov et"Zll2006l) . Thus, B oc R\j'^{lo/h)B^^ in Eq. ^ 
with Rm ^ 1, where Ru is the turbulent magnetic Reynolds num- 
ber of the outflow and Beq is the magnetic field strength corre- 
sponding to equipartition with the turbulent energy; Bcq — 5 ^lG 
in the Solar neighbourhood of the Milky Way. Thus, the outflow 
must be strong enough to support the dynamo action, Ru > 0.1 
according to our results (Figs [T] and |2ll. However, any outflow is 
removing the large-scale magnetic field as well, and thus adversely 
affects the dynamo action. Hence, the large-scale magnetic field de- 
cays when Ru > 0.8. The strength of the outflow optimal for the 
dynamo action is Ru ~ 0.3, and this is consisten t with the plau- 
sible range of U z in spiral galaxies estimated by Ishukurov et al.l 
(2006) asF^ = 1-2 km s^^ 

For Ru ~ 0.3, the dynamo achieves a steady state in about 
10^" yr for parameter values typical of the Solar neighbourhood 
of the Milky Way, with the amplification factor of about 10* in 
terms of magnetic energy. This imposes significant restrictions on 
the strength of th e seed magnetic field required for the dynamo 
( lBecketalJl993) . However, this estimate, often taken as represen- 
tative of spiral galaxies as a whole, only applies to the Solar vicinity 
of the Milky Way and, in addition, relies on various poorly known 
factors as well as on the approximations of this paper. Closer to the 
Galactic centre, the angular velocity of rotation fl oc r"^ is larger 
together with the dynamo number, D oc oc Q,^. For example, 
at a galactocentric distance equal to half the Solar orbit radius, the 
dynamo number is four times larger than near the Sun, and the dy- 
namo growth time r oc ("v/B — x/BT)"^ is more than five times 
shorter than near the Sun. This applies to other galaxies as well: 
for example, the growth rate of the large-scale magnetic field in the 
nearby galaxy M5 1 is estimated to be ten times larger than that in 
the Solar neighbourhood (§VII.9 in Ruzmaikin et al. 1988). 

These features of the mean-field dynamo facilitated by the ad- 
vective flux of magnetic helicity appear to be quite satisfactory. 
However, the steady-state strength of the large-scale magnetic field, 
about O.lBcq ~ 0.5 ^lG for D = 2Dc and Ru = 0.3 - see 
Eq. i20\ and Fig.|2l is n ear the low er end of the range observed in 
spiral galaxies, 1-5 /xG ( lBeckl2000l) . The main factor which makes 
the magnetic field strength low is B oc C~^^^ on lo/h 0.2 in 
Eq. J20b . The origin of this factor is the fact that the magnetic he- 
licity evolves over the time scale lo/rit — lo/vo - see Eq. {6} - 
which is shorter than the evolution time scale of the mean magnetic 
field, /rit- Thus, even in the presence of the advective helicity 
flux magnetic helicity can partially cancel the kinetic one before the 
large-scale magnetic field has grown enough. (Without any helic- 
ity flux, this cancellation is more complete and the mean magnetic 
field is catastrophically quenched.) The advection term in Eq. ^ 
opposes the growth of |am|, and so the steady-state strength of the 
mean field increases with Ru as (Ru/C)^^^. We note in this con- 
nection that it is important that Eq. (|5) involves a scale (/o) indepen- 
dent of the magnetic Reynolds number; otherwise, the strength of 
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the mean magnetic field would be catastrophically small for large 

m • 

Although the agreement with observations within an order 
of magnitude may be sufficient for a crude model explored, 
we discuss in the Sect. [5] an additional effect that can make 
the magnetic field stronger and accelerate its growth. We have 
considered an a d dition al flux of magnetic helicity, suggested by 
IVishniac & Choi bOOlh , which arises because of the symmetry 
breaking and anisotropy introduced in the turbulent flow by differ- 
ential rotation. This attractive mechanism is essentially nonlinear 
and can only be efficient if the large-scale magnetic field is strong 
enough - in fact, almost exactly as produced owing to the advective 
helicity flux (Sect. lSTt . Therefore, we believe that the primary role 
of the Vishniac-Cho helicity flux in galactic dynamos is to comple- 
ment the action of the advective flux of magnetic helicity. We have 
not explic itly included an other type of magnetic he l icity fl ux sug- 
gested bv lKleeorin et al.l 62000.) and lKleeorin et alj ( l2002h . which 
also relies on anisotropy of turbulence. However, the main elements 
of the functional form of that flux are, at least in the approximation 
used here, similar to that of the Vishniac-Cho flux, and we believe 
that the results would remain qualitatively unchanged. 

The additional helicity flux due to anisotropy of turbulence 
can be essential for supporting large-scale magnetic fields in star- 
burst galaxies or in galaxies with a relatively strong wind, where 
the strong outflow could otherwise suppress the dynamo action. 

Even for Ru = 0, the mean magnetic field achieves strength 
of order O.lBeq before it is catastrophically quenched. As fol- 
lows from Fig. |2l the field strength is maximum at a time of or- 
der 10^" yr, which is comparable to the galactic lifetime, and the 
subsequent decay by a factor of ten takes another lO^'' yr. This im- 
plies that some galaxies can have a significant mean magnetic field 
even if their fountain flow is too weak to alleviate the catastrophic 
quenching of the mean-field dynamo. 

A feature of both the earlier models of dynamical quenching 
and the present work is that the turbulent diffusion is assumed to be 
unaffected by magnetic field. Due to this reason, in the absence of 
a helicity flux, as a goes below a critical value (due to increase in 
I am I), turbulent diffusion leads to a decay of the mean field faster 
than at the resistive time-scale, resulting in the rapid decay of the 
mean field for Ru ~ mentioned in the previous paragraph. How- 
ever this is not accompanied in the models by an equally rapid de- 
cay of Qm. This is consistent with the fact that the steady-state so- 
lution of Eq. Qm ~ RynS ' B / B^q for Ru = 0, can be of 
order unity if f • B ~ R^^- One possibility is that current models 
are limited in how accurately they incorporate strict total helicity 
conservation. It is also possible that one may have a "turbulent" dif- 
fusive contribution to the small-scale helicity flux, when one goes 
to the next order in large-scale derivatives in Eq. [3] and the above 
problem might have resulted from the neglect of this flux. 

Alternatively, a mean-field decay without the correspondingly 
strong decay in Qm, can still be consistent with strict total helic- 
ity conservation if there is a preferential loss of the mean-field 
helicity through the boundary (without the corresponding loss of 
the small-scale helicity). Indeed such preferential loss of the large- 
scale field and its helicity (apparently due to the turbulent mag- 
ne tic diffusion of the mean field) are seen in the direct simulations 
of [Brandenburg & Doble3 ( 1200 ih where open boundary conditions 
were used. In these simulations helical turbulence was driven in a 
slab between open boundaries, with the field was taken to be purely 
vertical on the boundaries. The resulting a'^ -dynamo led to a mean- 
field, but with a steady-state magnetic energy decreasing with Rm. 
This presumably resulted from the fact that turbulent diffusion of 



the mean-field (and its helicity) through the boundary was more 
efficient than similar turbulent diffusive losses of the small-scale 
helicity. Such a behaviour is similar to the catastrophic quench- 
ing of the mean field that obtains in one-dimensional mean field 
dynamo models (including dyna mo models) with open bound- 
aries, such as this work or those of [Brandenburg & SubramairianI 
(2005c) and Shukurov et al. (200i), in the limit of zero small-scale 
helicity flux. Remarkably, in a domain with closed or periodic 
boundaries, where both large- and small- scale helicity losses are 
zero, a differe nt type of flm -depe ndent evolution obtains. Direct 
simulations of lBranden burg (2001) show that in this case the mag- 
nitude of the large-scale magnetic field can eventually reach larger, 
super-equipartition values, but the time scale to reach this value in- 
creases with Rm- Clearly, the case of open boundaries is more rel- 
evant for galaxies, and both the direct simulations and mean field 
dynamo models indicate that efficient removal of the small-scale 
magnetic helicity from the dynamo region is necessary for a healthy 
mean-field dynamo action. 

On the technical side, we have confirmed the applicability 
of the no-z approximation to the fairly complicated nonlinear dy- 
namo system discussed here. In particular, it is reassuring that, the 
time evolution of the magnetic field in the no-z model follows 
closely that obtained in numerical solutions of the corresponding 
partial differential equations, including the oscillatory behaviour 
( Shukurov et al. 2006). A nontrivial issue in this term approach is 
how to approximate J ■ B, since this vanishes in the no-z approxi- 
mation. Using a different approximation (Appendi x [At. we show 
that th is term is proportional to —[Dl^^^ BrB^. iKleeorin et al.l 
assume a similar form but without the factor ~\D\^^'^ which 
noticeably affects the qua ntitative results (especia lly because of the 
different sign of the term). lshukurov et alj ( l2006b approximate this 
term as being proportional to in their analytical estimate of 
the steady-state field, their Eq. (10), which affects the result for 
Ru ~ 0. With this form, the time evolution of the field (in the no-z 
approximation) is also not correctly reproduced for Ru = 0. 

The mechanisms of dynamo action discussed here can have 
extensive implications for galactic magnetic fields. For example, 
the modulation of the dynamo action by the outflow may contribute 
to the formation of magnetic arms, i.e., spiral-shaped regions of 
strong large-scale magnetic field located between the gaseous spi - 
ral arms, as in the galaxy NGC 6946 (e.g., Sect. 5 in lBeckll200a) . 
Magnetic arms occur not in all galaxies, and in some cases (e.g., 
M51) the large-scale magnetic field is stronger in the gaseous arms 
at some radii and between them elsewhere. It is reasonable to ex- 
pect that the intensity of the galactic fountain, quantified by Ru, 
is higher in the arms where star formation is more intense. If then 
the average value of Ru is less then 0.3, the enhancement of Ru 
in the gaseous arms will lead to a stronger magnetic field there. 
Otherwise, for stronger overall outflow, enhancement of Ru will 
suppress the magnetic field in the arms, which can produce mag- 
netic arms interlaced with the gaseous spiral arms. 



ACKNOWLEDGEMENTS 

We are grateful to E. Blackman, A. Brandenburg and J. Whitaker 
for useful comments. AS and KS were supported by the Lever- 
hulme Trust via grant F/00 125/N. SS would like to thank Council 
of Scientific and Industrial Research, India for financial support. 
AS gratefully acknowledges partial financial support of the Royal 
Astronomical Society. 



R 



© 2006 RAS, MNRAS 000,[T]l9] 



8 S. Sur, A. Shukurov and K. Subramanian 



APPENDIX A: PERTURBATION SOLUTION FOR THE 
atj-DYNAMO 

The kinematic aoj-dynamo in a thin disc is govern ed by the follow- 
ing e quations written in dimensionless form fe.g. jRuzmaikin et alJ 



ing_e 
[l988 



d , ^ . , d'^Br 

-Ra-^lyaBs) 
oz 



R^Br + 



dz^ 



with the vacuum boundary conditions at the disc surface, 

Brlz^l ~ Bd,\z = l — . 



(Al) 
(A2) 

(A3) 



where 7 is the growth rate of the magnetic field, dB/dt = 7B, 
and Bz can be recovered from the solenoidality condition. Here 
we derive an approximate solution of this eigenvalue problem for 
\D\ = \RaRui\ ^ 1. As it often happens with such asymptotic 
methods, the solution can be applied even for a relatively large \ D\. 
For the sake of definiteness, we assume D < 0, a(z) — sin nz and 
consider quadrupolar modes that are dominant in a thin disc. 

Using new variables = Ra^^r and = \D\~'^^^'B^, 
we rewrite Eqs iAU and iA2\ in terms of the dynamo number and 
then represented them in a symmetric form 



7B = {W + eV)B, 

where e — \D\^'"^ is a small parameter and for D < 0, 



(A4) 



W - 



VB 



oydz'^ 
d'^/dz^ 



o_ 

-Br 



-d{aB^)/dz 




and we have dropped dash at the newly introduced variables. 

The eigensolutions of the unperturbed (free-decay) system 
Anf>n = Wbn (wlth the above boundary conditions) are doubly 
degenerate and given by 
2 



An 



bri = 



b' — 



-\/2 cos[7r(n + ^)z] 




V2cos[n{n + ^)z] 



(A5) 



where we have normalized them to b^ dz — b'^ dz — 1; the 
eigenfunctions should not be confused with the small-scale mag- 
netic field denoted b in the main text. 
The expansions 



7 = 70 + £71 + <s 72 + • ■ • , 

B = Cobo + C'ob'o + eC-ibi + ^c[b[ + . 



(A6) 
(A7) 



are substituted into Eqs dAU and iA2\ . terms of like order in e col- 
lected, the dot product of the resulting equations taken first with 6„ 
and then with b'„, and results are integrated over ^ z ^ 1. To the 
lowest order, this yields 70 = Aq. A homogeneous system of alge- 
braic equations for Co and Cq follows from terms of order e, whose 
solvability condition yields 71 = (Vo'oVqo')^'''^ ~ iV^' 
Co = C071/V00' = -2Co/v^, where Vnm = bn ■ Vbm dz 
are the perturbation matrix elements, whose direct calculation 
yields Voo = Vo'c = Vio = Vvo = Vj/o' = 0, Vqo' = -tv/A, 




Figure Al. The approximate eigenfunctions Br (dashed) and (solid) 
from Eqs (A8) and <A9) for D = -8, -20, -30 (with Ki -8), Cq = 1 
and Ra = 1. 



^o'o = ~1 ™d Vio' = — 37r/4. To the same order in e, we sim- 
ilarly obtain inhomogeneous algebraic equations for Ci and Ci, 
which yield Ci = and Ci = C7oVio'/(Ao-Ai) = 3Co/(47r='/2). 

In terms of the original variables Br,^, the resulting approxi- 
mate solution is given by 



Br 



Bi 



RaCo 



-2Co\ 



(A8) 



(A9) 



(AlO) 



where Co is an arbitrary constant. Thus, 7 > for _D < Dc 
with Dc ~ — 7r^/4 ~ —8. This estimate of is imp ressively 
close to that obtained numerically jRuzmaikin et al.ll988[) . We plot 
the components of the magnetic field for various values of the dy- 
namo number in Fig. lAll the result reproduces very closely the 
more accurat e nurnerical solutions presented, e.g., in Fig. VII. 1 of 
iRuzmaikin et a lJ (l988l) . including such a subtle detail as a zero of 
Br at < z < 1 for D < Dc which shifts to smaller 2 as |D — Dc] 
increases. 

The current helicity density of this magnetic field follows as 



J-B^^^^\D\ 



'^'BrB, 



(All) 



to the lowest order in D. This estimate is used in Sect.|3] 
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